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Abstract 



The partition function and specific heat of a system consisting of a finite 
number of bosons confined in an external potential are calculated in differ- 



> 

in 
0^ 
O 

Tij- . ent spatial dimensions. Using the grand partition function as the generat 

o , 

. ing function of the partition function, an iterative scheme is established for 

"j^ . the calculation of the partition function of system with an arbitrary number 

I . of particles. The scheme is applied to finite number of bosons confined in 

^ , isotropic and anisotropic parabolic traps and in rigid boxes. The specific heat 

O ■ 

as a function of temperature is studied in detail for different number of par- 
^ . tides, different degrees of anisotropy, and different spatial dimensions. The 

cusp in the specific heat is taken as an indication of Bose-Einstein condensa- 
tion (BEC). It is found that the results corresponding to a large number of 
particles are approached quite rapidly as the number of bosons in the system 
increases. For large number of particles, results obtained within our iterative 
scheme are consistent with those of the semiclassical theory of BEC in an 
external potential based on the grand canonical treatment. 
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I. INTRODUCTION 

Following the foundation laid down by EinsteinS, discussions on the phenomena of Bose- 
Einstein Condensation (BEC) have been of constant interest to physicists and can be found 
in standard textbooks on statistical physics.H The conventional treatment is usually based on 
the grand canonical ensemble within the context of an ideal Bose gas and the thermodynamic 
limit. The recent beautiful experimentsi on BEC in atomic gases, however, lead to the need 
of a deeper understanding of the problem of BEC in confined systmes with finite number of 
bosons. The fact that the number of bosons is finite and relatively small in the trap created 
experimentally and that the ultracold atomic gas does not exchange particles with reserviors 
implies that the canonical ensemble is more appropriate.Bi 

It is interesting to study how the thermodynamical quantities such as the specific heat 
evolves as the number of particles in the system increases. In Fermion systems, the differ- 
ences among results obtained within the different ensembles have been studied.§ Recently, 
there have been some work on finite boson systems in a parabolic confinement within the 
context of micro-canonical and canonical ensembles.0i In the present work, using the grand 
partition function as the generating function of the partition function, an iterative scheme is 
established for the calculation of the partition function of system with an arbitrary number 
of particles. The scheme is applied to bosons confined in isotropic and anisotropic traps 
and in rigid boxes in arbitrary spatial dimensions. In particular, we aim out studying the 
effects of dimensionality, anisotropy, and finite number of particles. Although there exist 
many interesting approaches to the problem in the literature,! we concentrate on comparing 
our results, wherever possible, with those in the conventional theory and in the semiclassical 
theory of Bagnato et a/.,i which is basically an extension of the conventional theory to the 
case with confining potentials. 

The plan of the paper is as follows. In Sec. II, the general formalism and the recursive 
scheme are given. Results for bosons confined in isotropic and anisotropic parabolic po- 
tentials are given in Sec. III. The corresponding results for confinement in rigid boxes are 
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given in Sec. IV together with a discussion on the vahdity of the integral approximation 
commonly used. A brief summary is given in Sec. V. 



II. FORMALISM 

We consider a system of non-interacting bosons. The partition function Z]sr(v,T) is 
related to the A^-th derivative with respect to the fugacity z of the grand partition function 
Q(^z,v,T) through 

1 d''Q{z,v,T) 
^^(^' ^) = ]v! dl^ 1^=°' 

where T is the temperature and z = exp{^/kT) is the fugacity. Here v is the variable 

representing the volume of the system for confinement in a rigid box and the confining 

strength in the case of a parabolic confining potential. In general, the grand partition 

function for an ideal bose gas can be written aM 

where the product is over all the single particle states labelled by p with energies ep, and 
f3 = 1/kT with k being the Boltzmann constant. The zeros of the chemical potential and 
the single particle energy are set to be the same. Combining Eqs.(l) and (2) gives a relation 

1 ^ 

= ]^ XI Bj^N-j (3) 
i=i 

between the partition function Zn of a A^-particle system and the partition functions Z^/, 
with A^' = 1,---,A^ — 1, for systems with less than A^ particles. The factor Bj is defined as 
Bj = J2p^~''^^''- For simplicity, the arguments v and T are left out in the partition functions 
in Eq.(3). The internal energy E and the heat capacity Cy can then be calculated using the 
standard formulae E = —din Z]\f/dj3 and = dE/dT. These expressions involve the first 
and second derivatives of the partition function with respect to temperature. 

The first and second derivatives of the partition function of a A^-particle system can also 
be related to quantities in systems with less particles. Let 6 be a dimensionless parameter 
representing the temperature. From Eq.(3), we have 
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dZN-j dBj 



(4) 



and 



d'^Z 



N 



N 



N 

E 



'd-'ZN-j^^dBjdZM-j^d-'B 



(5) 



db^ db db db^ 

Equations (3)-(5) imply that the partition function and other thermodynamical quantities 
of a A^-particle system can be calculated from the partition functions and their derivatives 
of systems with less particles iteratively by using the conditions Zq — 1, dZ^/db — 0, 
and d^Zo/db"^ = 0. These conditions, when substituted into Eqs.(3)-(5), give the correct 
single particle partition function Zi and its derivatives. It should be pointed out that these 
conditions are independent of the form of the confining potential and spatial dimension. 
The effects of the confinement are reflected in the energy spectrum and the explicit form of 
Eqs.(3)-(5) when they are applied to specific problems. In what follows, we apply Eqs.(3)- 
(5) to systems with finite number of non-interacting bosons confined either in a parabolic 
potential or in a rigid box. In particular, the position of the cusp in the heat capacity as 
a function of temperature is taken to be an indication of Bose-Einstein condensation in the 
system. We study the dependence of the position of the cusp on the number of particles in 
the system. 



III. PARABOLIC CONFINEMENT 

A. Isotropic confinement in 1, 2, and 3 dimensions 

For M-dimensional isotropic parabolic confinements characterized by the energy huj, the 
energies of the single particle states are given by 

M 
i=l 

where rii {i = 1, ■ • ■ , M) takes on non-negative integers and the zero-point energy has been 
absorbed in the definition of the zero of energy. Identifying the dimensionless parameter b 
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as 6 = exp(— 1/r), where r = kT/hu is the effective temperature, and using Eqs.(3)-(5), we 
have, in M dimensions, the relations 
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db 



N-j 



S (1 - b^) 



M ■ 
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1 ^ 1 
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., j-l + (jM + l)6^- 

(1 - b^y 



1 - b^ 
Zn-j 



db 



(9) 



Equations (7)-(9), together with the conditions on Zq and its derivatives, can be used itera- 
tively to obtain Zi, ■ ■ -, Z^ and their derivatives. These equations can be easily implemented 
numerically for any spatial dimensions. The heat capacity can be evaluated using 
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Zn db \ Zn db J Z^ db"^ 

It should be pointed out that, although exact expression for Z^ may not be possible for 
general dimensions, exact expressions exist for Z^ and in one dimension (ID) and they 
read 



7(m 
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nti(i 



(11) 



and 



'huVSL fb' 



kT) ,5(1-6^)2- 

Slightly different expression for Z^^^^ has been obtained by Brosens et a/..0 The difference 
comes from our choice of setting the zero of energy at the ground state energy of a m- 
dimensional harmonic oscillator. This choice, of course, does not affect the result of the heat 
capacity. Figure 1 shows the results of the specific heat C^/Nk calculated using Eq.(lO) as 
a function of kT / {TlujN^/^'') for systems in M = 1, 2, and 3 dimensions (ID, 2D, 3D) with 
different number of bosons. In Fig. 1(a), it is obvious that the results obtained numerically 
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using Eqs.(7)-(10) agree well with the results obtained from the exact expression given by 
Eq.(12). This, in turn, indicates that our iterative approach is reliable. In ID, the specific 
heat approaches its large limit very rapidly in that the curves corresponding to systems 
with = 10 and = 10^ particles are very close to each other. As temperature increases, 

approaches its limiting value of Nk. In contrast to 2D and 3D systems, varies smoothly 
with temperature and there is no cusp. This implies that there is no BEC in ID system with 
parabolic confinement. This result is consistent with that of the semi classical treatment of 
Bagnato and Kleepner.il These authors studied BEC in traps within the grand canonical 
ensemble and the thermodynamic limit so that the notion of a density of states is valid. 
They found that ID ideal Bose gas will display BEC only if the external potential is more 
confining than the parabolic form. 

Figure 1(b) shows that cusp appears in the specific heat when the number of particles in 
the system is large enough in 2D isotropic parabolic confinement. The position of the cusp 
shifts to higher temperature as the number of particles increases. As there will be, strictly 
speaking, no sharp transition in a system with finite number of particles, the appearance 
of a cusp in the specific heat is taken to be an indication of BEC. For system with A^ = 
10^ particles, the cusp is at a temperature kT/ {htoN^^'^) ~ 0.75; while the semiclassical 
treatment]! gives a transition temperature Tc at kTc/ (hcuN^^'^) = \^/n ^ 0.7797 for A^ oo. 
Thus, our results are consistent with the semiclassical results and we expect the position of 
the cusp in Cy to approach a value slightly larger than kT / {TiujN^^'^) ~ 0.75 as A^ increases. 

In 3D isotropic parabolic confinement, the specific heat shows a cusp as a function 
of temperature for systems with large enough A^ with the position of the cusp shifts up in 
temperature as A^ increases and approaches a limit (see Fig. 1(c)). However, the behaviour is 
quite different from that in 2D. As A^ increases, the cusp becomes sharper and decreases 
rapidly across the cusp and eventually leads to a discontinuity in C^. For systems with 
A^ = 10^ particles, the position of the cusp is at kT/{hiuN^^^) 0.9 and the discontinuity 
in is about — 6.5A^A;. According to the semiclassical treatment for BEC in 3D harmonic 
traps,i the transition temperature is at kTc/ {htuN^^^) ~ 0.912 for A^ oo, and there will 
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be a discontinuity in of the amount given by [Cv{T^) — Cv{T~)]/Nk —6.57. Therefore, 
our results for large N are consistent with that of semiclassical theory. In contrast, the 
specific heat in is continuous across the cusp in 2D. 

B. Anisotropic confinement in 2 dimensions 

We choose to illustrate the effects of anisotropy by considering two-dimensional systems. 
The energies of the single particle states in a 2D anisotropic harmonic confinement are 

£ni,n2 = + «"'2), 711,^2 = 0, 1,2, •• •, (13) 

where a — 002/^^1 is the ratio of the frequencies of the harmonic potentials in the two 
directions. In this case, the factor Bj is defined to be Bj{b) = 1/(1 — 6^)(1 — b'^^), where 
b = exp(l/T) with r = kT/hcui. Equations (3)-(5) can then be applied iteratively to obtain 
the heat capacity. Figure 2 shows the dependence of the specific heat on temperature by 
plotting Cy/Nk as a function of T/{aNy/'^ = kT /n{ijJiUJ2Y''^N'^''^ for different values of a 
with different number of particles in the system. Figure 2(a) shows the case with fixed 
degree of anisotropy a = 80. For = 10^ particles in such an anisotropic trap, the 
cusp appears at t/(q;A^)^/^ ^ 0.71; while semiclassical theory gives a value of 0.7797 for 
N ^ 00. For small number of particles, say N = 10, Cy appears to approach the ID 
limit of Nk at low temperature. However, as temperature increases, increases again and 
eventually approaches the 2D hmit of 2Nk at high temperature. This behaviour can be 
understood in that the anisotropy leads to different energy spectra in the two directions. At 
low temperature, only the levels corresponding to the less confining direction are populated 
and thus the system appears to behave in a one-dimensional way. At higher temperature, 
the levels corresponding to the larger confining direction start to become populated and the 
system takes on two-dimensional character. Figure 2(b) shows the effects of anisotropy for 
a fixed number of particles N — 10^ in the system. For small a, there is no cusp in - a 
character of a ID system. However, it should be pointed out that even for small values of 
a, Cy approaches the 2D limit of 2Nk at high temperatures. 
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IV. RIGID BOX 



The iterative approach can also be used to calculate the thermodynamical quantities 
of a system of N non-interacting Bosons confined in a M-dimensional rigid box of volume 
V = . The energies of the single-particle states for a particle of mass m in such a box are 
given by 



(14) 



where the quantum numbers rii {i = 1,2, ■■■ , M) takes on positive integers. The factor Bj in 
M dimensions is given by Bj{h) = [Qj{b)]^^ with Qj{b) being the factor Bj in ID defined by 
0j(^) = J2'^=i b^"'^ ■ The dimensionless parameter b is taken to be 6 = exp(— 1/r), where r is 
the effective temperature given by r = SmkTL'^/h'^. In m dimensions, Eqs.(3)-(5) become 

1 N 

(15) 
(16) 
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^3 



0f db 

1 d^Q, 
^ Wj~dh^ 



db 



(17) 



The function Oj(&), which is related to the Jacobi's Theta-function,i and its derivatives 
can be evaluated numerically by summing the corresponding infinite series to convergence. 
The heat capacity can then be obtained using Eq.(lO). Figure 3 shows the specific heat for 
different number of bosons in a rigid box in ID, 2D, and 3D. Figure 3(a) shows that in ID 
rigid box the results for systems with a small number of particles approach that of a larger 
number of particles quite rapidly. This feature is similar to that in a ID parabolic trap. 
Another feature is that the high temperature limit of = 0.5k is approached quite slowly. 
For example, C^/Nk increases from approximately 0.478 to 0.489 as t/N"^ increases from 20 
to 80. 
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Figure 3(b) shows that in 2D the large behaviour is approached very rapidly. The high 
temperature limit is approached more rapidly than in the ID case. There is no cusp in the 
specific heat for both ID and 2D. This is consistent with the well-known result that BEG 
for ideal bose gas in the absence of external confining potential does not exist in dimensions 
lower than three. Figure 3(c) gives the results in 3D. The position of the cusp shifts to 
approach the large N limit from above as the number of particles increases. For = 10'^, 
our numerical results give the maximum specific heat at the cusp the value C^/Nk = 2.087 
at the corresponding temperature SmkTL'^ /h'^N'^^^ ^ 0.7556. It should be contrasted with 
the standard ideal Bose gas results in 3D in the thermodynamic limit that the maximum 
specific heat is 1.925A^A; at the transition temperature given by SmkT^L^/h'^N^/^ ^ 0.6713. 
In contrast to the results corresponding to a 3D parabolic potential, the specific heat is 
continuous across the cusp. 

The infinite series in Qj{b) is sometimes approximated by an integral! which yields Qj ~ 
(7rr/4j)^/^. With this approximation, a suitable dimensionless parameter is 6 = (7rr/4)^^/^. 
An approximate expression of the the factor Bj{b) in M dimensions is 

B,{b)^r''^%. (18) 

IFtoiu. the identity 

j2 e-^"'/- = -{—y/^ + (—y/^ e-^"""^/^ - -, (19) 

n=l 2 j j g^-^^ 2 

the integral approximation amounts to retaining only the first term. Within this approxi- 
mation, Eqs.(15)-(17) read 

1 ^ 

ZN = j:j:r''^'bZ^-^^ (20) 



dh Nj^{ 

q2 7 1 TV 



(21) 
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(22) 
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Within this approximation, the formula for heat capacity is 
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and the calculated results are shown in Fig.4. The general features are quite similar to that 
of the exact results shown in Fig. 3. It can be shown that the integral approximation gives 
reliable results except for systems with only a few particles, in those cases the approximation 
may lead to erroneous results. 



V. CONCLUSIONS 

We have studied the BEC of finite number of confined bosons. An iterative scheme has 
been used to obtain the specific heat for systems with different number of particles confined 
in an external potential. As emphasized by van Hove,0 no singularities would appear in 
the partition function and the thermodynamical quantities of a finite system. Therefore 
the cusp of speific heat curve was taken as an indication of BEC. Numerical results show 
that the specific heat approaches its large N limit rapidly as the number of particles in 
the system increases. Results in different spatial dimensions are obtained. Results for 
systems with large A^, typically A^ ~ 10'' in our calculations, are consistent with those in 
the semiclassical theory of Bagnato et a/..i As the early experiments have typically 10^-10^ 
particles in the trap,! we expect that treatments within the canonical and grand canonical 
ensembles would give at least qualitatively consistent results. Although the iterative scheme 
is applied here to Bose systems, it can equally be applied to systems with finite number of 
Fermions confined in an external potential. In this case, the variable Bj should be modified 
to be Bj = Ep{-iy~^e-^f^'^. 
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FIG. 1. The specific heat of a finite number N of bosons confined in an isotropic parabohc 
potential as a function of the dimensionless temperature rjN^I^ = kT/{TiLoN^/^), where M is 
the spatial dimension for (a) M = 1, (b) M = 2, and (c) M = 3. The lines are results of the 
iterative scheme for different values of N. The symbols in (a) are the results obtained using the 
exact expression Eq.(12). 



14 



4.0 



2.0 



- N=10 
- -N=100 

9 ^ - - N=1000 

I ^'h ^--^ N=10000 

P / . . -A -A -A^ -A^ -A-A 



,A' 



A 



i-^ 0=80 



_ / / 

7 ^''^ 

A/ ^ 



(a) 



0.0 ^ 
4.0 r 



a=10"' 
a=10"^ 

/ \ ^ a=10" 



^ 



2.0 



'to 



N=10000 



4) 

/ (b) 



4 

0.0 



0.0 1.0 2.0 3.0 

1/2 



x/(aN) 



15 



FIG. 2. The specific heat of N bosons confined in a two-dimensional anisotropic parabohc 
potential as a function of the dimensionless temperature T/{aN)^/^ = kT /%{uiU2Y^'^N^^'^ ■ (a) 
The degree of anisotropy a is taken to be 80 and the number of particles takes on four different 
values oi N = 10, 10^, 10^, and 10^ . (b) The number of particle is fixed at = 10^, and the 
degree of anisotropy a takes on a = 10~^, 10~^, and 10~^. 
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FIG. 3. The specific heat of bosons confined in M-dimensional rigid box as a function of 
the dimensionless temperature r/iV^/^ = SmkTL'^ /h?N'^/^ for (a) M = 1, (b) M = 2, and (c) 
M = 3. Results are obtained using Eqs.(15)-(17) by summing up the infinite series for @j{b) to 
convergence. 
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FIG. 4. The specific heat of bosons confined in M-dimensional rigid box as a function of 
the dimensionless temperature t/N'^/^ = SmkTL'^ /h?N'^/^ for (a) M = 1, (b) M = 2, and (c) 
M = 3. Results are obtained using Eqs.(20)-(22) within the integral approximation. 



20 



